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ABSTRACT 

In this Letter, we present a systematic study of lenticular (SO) galaxies based on mid-infrared 
imaging data on 185 objects taken using the Spitzer Infra Red Array Camera. We identify the SOs 
hosting pseudobulges based on the position of the bulge on the Kormendy diagram and the Sersic 
index of the bulge. We find that pseudobulges preferentially occur in the fainter luminosity class 
(defined as having total if-band absolute magnitude Mk fainter than —22.66 in the AB system). 
We present relations between bulge and disk parameters obtained as a function of the bulge type. 
The disks in the pseudobulge hosting galaxies are found to have distinct trends on the r c — rd and 
A*d(0) — r d correlations compared to those in galaxies with classical bulges. We show that the disks 
of pseudobulge hosts possess on average a smaller scale length and have a fainter central surface 
brightness than their counterparts occurring in classical bulge hosting galaxies. The differences found 
for discs in pseudobulge and classical bulge hosting galaxies may be a consequence of the different 
processes creating the central mass concentrations. 

Subject headings: galaxies: photometry — galaxies: formation — galaxies: fundamental parameters 



1. INTRODUCTION 

Bulges are a central piece of the puzzle of galaxy for- 
mation and evolution. Recent work - theoretical and 
observational - has revealed that bulges come in two fla- 
vors: those thought to have formed through violent pro- 
cesses such as hierarchical clustering or major mergers 
and referred to as "classical bulges" and those thought 
to have formed through secular processes and referred to 
as "pseudobulges" (Kormendy & Kennicutt 2004). Many 
studies have been carried out proposing and examining 
various criteria for identifying these two classes and dis- 
tinguishing between them. For example, Carollo et al. 
(1998) identify pseudobulges based on the presence of nu- 
clear structure as seen in images taken using the Hubble 
Space Telescope (HST). Fisher k Drory (2008) employed 
a similar method and proposed a classification of bulges 
into "classical" and "pseudo" based on the Sersic index 
of the bulge. Carollo et al. (2001) find a difference in 
the average optical-near-infrared V — H color of the two 
kinds of bulges. The two types exhibit different behavior 
on various well known correlations between the struc- 
tural parameters of the galaxy. The studies also show 
a smooth (rather than sharp) transition from one type 
of bulge to the other, possibly implying the existence of 
bulges with a mixture of properties (Gadotti 2009; Fisher 
& Drory 2010). 

The above studies have focussed on bulge dichotomy 
using samples comprising of all galaxy types with a bulge 
and a disk but it can also be interesting to study this di- 
chotomy within the same morphological class. In the 
present work, we consider SO galaxies which form an in- 
termediate transition class between ellipticals and spi- 
ral galaxies on the Hubble tuning fork diagram (Hub- 
ble 1936). Recent studies hint at existence of sub- 
populations within SO galaxies, with different properties 
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and formation histories. Barway et al. (2007, 2009, 2011) 
have proposed a division of SO galaxies into bright and 
faint classes using a dividing if-band absolute magnitude 
of —24.5 in the Vega system. The properties and behav- 
ior with respect to various well known correlations be- 
tween morphological parameters are different for the two 
classes. For example, Barway et al. (2007) find a strong 
positive correlation between the bulge effective radius r c 
and the disk scale length for fainter galaxies while a 
weak anti-correlation is found in the case of bright galax- 
ies. Barway et al. (2009) report distinct trends for bright 
and faint SOs on correlations such as the Kormendy rela- 
tion, the photometric plane etc. The differences can be 
interpreted as the dependence of the formation history 
of the galaxies on the luminosity and the environment of 
the galaxies: the more luminous SO galaxies should have 
formed through more violent processes while the faint 
ones through secular evolution. 

If pseudobulges in SO galaxies have formed through 
processes wherein the disk material rearranges itself to 
form the bulge, it is expected that the disk itself un- 
dergoes a change in its properties when giving rise to 
the bulge. In this Letter, we have carried out a mid- 
infrared study of a large sample of SO galaxies to exam- 
ine the properties of the bulges and the host disks. The 
Letter is organized as follows. Section 2 describes the 
sample selection and the data. Section 3 describes the 
various results found and Section 4 summarizes the find- 
ings and their implications. Throughout this Letter, we 
have assumed the standard concordance cosmology i.e. 
H Q = 70 kms" 1 Mpc" 1 , ti m = 0.3 and fi A = 0.7. 

2. SAMPLE AND DATA ANALYSIS 

To assemble a statistically significant sample of SO 
galaxies, we started with the RC3 catalog (de Vau- 
couleurs et al. 1991) and identified all galaxies classi- 
fied as SO (with numerical Hubble stage —3 < T < 0). 
The RC3 contains 3657 such galaxies. We impose a mag- 
nitude cut with total apparent B magnitude Bt < 14.0 
which reduces the sample size to 1031 galaxies. For these 
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galaxies, we obtained the total apparent K s magnitudes 
from the Two Micron All Sky Survey (2MASS) All Ex- 
tended Source Catalog and converted these to absolute 
magnitudes using redshifts obtained from the NASA Ex- 
tragalactic Database (NED)0 

The mass content of a galaxy is dominated by low mass 
stars and radiation at infrared wavelengths takes into ac- 
count their contribution and is also relatively free from ef- 
fects of dust extinction. It is therefore appropriate to use 
data at near- or mid-infrared wavelengths for the study 
of galaxy morphology. We used the imaging data at 3.6 
/xm taken using the Infrared Array Camera (IRAC) on 
board the Spitzer Space Telescope. IRAC offers deep im- 
ages at mid-infrared wavelengths free from the problems 
of sky emission that plague ground based observations. 

We crossmatched the above magnitude limited sample 
with the data available in the Spitzer Heritage Archive 
(SHA0 and found that imaging data for 247 galaxies 
were available. Spitzer takes dithered observations of the 
extended objects, which were downloaded from the SHA 
and coadded to a final mosaic using the Mosaicking and 
Point Extraction (MOPEX) tool provided by the Spitzer 
Science Centre. MOPEX performs the necessary prepro- 
cessing steps which include accounting for optical dis- 
tortions, image projection and outlier pixel rejection to 
create a final science mosaic. 

In order to derive the structural parameters of vari- 
ous galaxy components, we employed the technique of 
two-dimensional decomposition of galaxy light. We used 
the program GALFIT (Peng et al. 2002) for this purpose. 
In the first run, we fitted all galaxies with a bulge and 
a disk component using a Sersic (Sersic 1968) and an 
exponential profile respectively. For galaxies where the 
residual image obtained by subtracting the point spread 
function convolved best-fit model from the observed im- 
age revealed a bar, a second run of fitting was performed 
by adding another Sersic component to describe the bar. 
The bulge properties thus determined are free from the 
systematics that arise if a bar which is present remains 
unaccounted (Gadotti 2008; Laurikainen et al. 2005). 
During the course of the decomposition, galaxy images 
with poor signal-to-noise ratio and/or complicated mor- 
phological features such as tidal tails due to recent merg- 
ers were discarded from the sample. The final sample 
comprises 185 galaxies with a median redshift of ~ 0.005 
and a standard deviation of ~ 0.002. Further, the sam- 
ple size is limited only by the availability of good qual- 
ity data for non-interacting / undisturbed SOs in the 
SHA. The SOs were classified as bright and faint fol- 
lowing the criterion by Barway etal. (2007, 2009, 2011), 
who proposed a dividing Vega magnitude of —24.5 on 
the 2MASS K s absolute magnitude scale for this pur- 
pose. We transformed this division line to the AB sys- 
tem using i^s(AB) = K s (Vega) + 1.84 as suggested by 
Muhoz-Mateos etal. (2009). The sample comprises 37 
bright and 148 faint galaxies as determined using the 
new dividing magnitude of —22.66. 

We have used the AB magnitude system (Oke & Gunn 
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Fig. 1. — Kormendy relation with filled and open circles rep- 
resenting bright and faint SOs respectively. The solid line is the 
best-fit line to Coma cluster ellipticals while the dashed lines mark 
the 3-ct limits. 

1983) throughout this Letter and where necessary, mag- 
nitudes obtained from various catalogs have been trans- 
formed to this system. 

3. PSEUDOBULGES IN SO GALAXIES 
3.1. Identifying Pseudobulges 

Fisher & Drory (2008) classified as pseudobulges those 
bulges which had the presence of nuclear spiral arms, 
bars and/or rings as seen in the high resolution near 
V - band images taken using the Hubble Space Telescope 
(see also Carollo et al. 1998). The rest were classified as 
classical bulges. They found that classical bulges have 
Sersic index n > 2 and follow the same relations be- 
tween various structural parameters as elliptical galaxies. 
But these relations were not obeyed by the pseudobulges 
which are found to have n < 2. 

It appears from above that the Sersic index can be used 
to classify the bulges into the two types, classifying the 
bulges with n < 2 as pseudobulges and those with n > 2 
as classical. If we use this criterion to classify the bulges 
of the SO galaxies in our sample into these two types, 
we have 111 classical bulges and 74 pseudobulges. Only 
a small fraction (13/74) of pseudobulges belong to the 
bright SOs with the rest (61/74) being associated with 
faint galaxies. However, the typical error bars on n can 
be quite significant. Our experience with various data 
shows that this error can be as large as 20% while Gadotti 
(2008) reports errors as large as 0.5. The parameters of 
the Sersic profile n and r e are coupled which leads to 
an additional uncertainty in n (see Trujillo et al. 2001 
and references therein) . Thus classifying bulges using n 
alone can lead to ambiguity. 

Following Fisher & Drory (2008), we tried to use high 
resolution HST images for constraining the classification 
of the bulges by searching for evidence of nuclear struc- 
ture. For our SO sample, high resolution HST imaging 
data were available for only 110 of the 185 galaxies in 
various optical bands. In the absence of a homogeneous 
sample of HST images, it was not possible to impose 
constraints using this method. 

A different approach used by Gadotti (2009), which 
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avoids the difficulties associated with using n alone, in- 
volves classifying bulges based on their position on the 
Kormendy diagram (Kormendy 1977). This is a plot 
of the average surface brightness of the bulge within its 
effective radius (/ib(< r e )) against the logarithm of the 
effective radius r c . The elliptical galaxies are known to 
obey a tight, linear correlation on this diagram. Classi- 
cal bulges resemble elliptical galaxies and obey a similar 
relation while pseudobulges lie away from it. Any bulge 
which deviates by more than three times the rms (root 
mean square) scatter from the best-fit relation for ellip- 
ticals is classified as a pseudobulge by Gadotti (2009). 

The Kormendy plot for our sample is shown in Fig- 
ure [1] with filled and empty circles respectively denoting 
bright and faint galaxies. Since the classification of a 
bulge depends on its deviation from a best-fit line to the 
ellipticals, the decomposition data from a study by Khos- 
roshahi et al. (2000), of Coma cluster ellipticals done in 
the if-band, was used for comparison with our data. The 
magnitudes were transformed to the 3.6 /im band using 
the relation K — m^^ = 0.1 (Toloba et al. 2012; Falcon- 
Barroso et al. 2011) and to the AB magnitude system by 
adding 1.84 following Munoz-Mateos et al. (2009). The 
best-fit line for these ellipticals is 

( /ib (< r e )) = (2.567 ± 0.511) log r c + 16.595 ± 0.296. 

The rms scatter in (/J.b(< r c )} for the Coma cluster ellipti- 
cals is 0.592. The best-fit line is shown in Figure [TJ The 
dashed lines in the diagram enclose a region where all 
points within three times the rms scatter can be found. 
Following Gadotti (2009), points below this region can be 
classified as pseudobulges and there are 50 such points. 

Interestingly, most pseudobulge hosting SOs thus clas- 
sified are faint. Only three of the 37 bright galaxies are 
classified as pseudobulge hosts based on their position on 
the Kormendy diagram. If we assume that both bright 
and faint SOs have the same fraction (32%) of pseudob- 
ulges, we can estimate the probability of finding only 
three pseudobulges among bright SOs using the binomial 
distribution. This probability is found to be ~ 10~ 4 , 
suggesting that it is highly unlikely that bright and faint 
galaxies have the same chance of hosting a pseudobulge. 

A set of prescriptions has been suggested by Kormendy 
& Kennicutt (2004) to put constraints on the bulge classi- 
fication. The more constraints a bulge satisfies, the more 
secure its classification becomes. Apart from examining 
HST images for various nuclear structures, one could in- 
vestigate the D n (4000) indices (Gadotti 2009) to check 
for recent star formation as pseudobulges are known to 
contain a recently formed stellar population and may 
even be forming stars actively. However, Z? n (4000) mea- 
surements, which are obtainable from the Sloan Digital 
Sky Survey (SDSS) data products, are again not avail- 
able for a large fraction of the SO galaxies in our sample. 
The SDSS is constrained to a specific portion of the sky 
while our sample is an all-sky sample. 

In the absence of the above indicators to have more 
secure classification of bulges in our SO sample, we chose 
to impose a cut off on the Sersic index n for the pseudob- 
ulge candidates obtained from the Kormendy diagram as 
explained above. We accept as a pseudobulge only those 
bulges which lie below the 3<r line and have n < 2. This 
reduces the final sample of pseudobulges to 27 (of which 
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Fig. 2. — Plot of r vs with filled and empty circles denoting 
classical bulges (CB) and pseudobulges (PB), respectively. 
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Fig. 3. — Scaled histogram showing the distribution of logrj. 
The solid and dotted lines represent distributions for classical and 
pseudobulges, respectively. 

two are bright) . All other bulges are classified as "classi- 
cal" . Where available, we have used HST images to con- 
firm the presence of nuclear structure in these galaxies. 
Further, wherever possible, we have checked HST images 
of those galaxies classified as classical bulge hosts and 
verified that these galaxies do not show any nuclear fea- 
tures. We believe that our pseudobulge identification is 
reasonably secure. The subsequent sections are based on 
this sample of 27 pseudobulges and 158 classical bulges. 

3.2. Disk Correlations as Function of Bulge Type 

If pseudobulges are formed by processes within the 
disk, one expects a correlation between the bulge and 
disk scale lengths (Courteau et al. 1996), which is pre- 
dicted by secular bulge formation models (Martinet 1995; 
Combes 2000; see Carollo et al. 1999 for comprehen- 
sive review articles). Also, one expects an imprint of the 
pseudobulge formation processes on correlations involv- 
ing disk parameters. 

In Figure [5] we plot bulge effective radius r c as a func- 
tion of the disk scale length r^- The filled and empty cir- 
cles respectively denote classical bulge hosting SOs and 
pseudobulge hosts. The empty circles denoting the pseu- 
dobulge hosts are shifted toward the left relative to the 
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Fig. 4. — Plot of disk central surface brightness as a function its 
scale length. The filled and empty circles represent the classical 
and pseudobulges, respectively. The best-fit straight lines to them 
are the solid and the dashed lines. 



filled circles which represent classical bulges. This in- 
dicates that the disk scale length is smaller on average 
for pseudobulge hosts than classical bulge hosts. This 
can be seen in Figure [3] The mean disk scale length 
for pseudobulge hosts is 1.6 kpc while that of classical 
bulges is 4 kpc. A t-test rules out at greater than 99.9% 
confidence, the possibility of this difference arising from 
random chance alone. This indicates that pseudobulges 
preferentially occur in smaller disks. The nature of the 
bulge depends on the formation history, in particular on 
whether the bulge is a product of a process such as a 
major merger or whether it has formed from secular evo- 
lution. If we assume that there is no reason to believe 
that only galaxies with larger disks should undergo such 
a merger it might be possible that the lower scale length 
of the disk in the case of a pseudobulge is due to the 
processes within it that grew the bulge. 

If secular processes that grow pseudobulges leave be- 
hind an imprint on the progenitor disks, one should be 
able to find them on correlations involving disk parame- 
ters. In Figure HI we plot the central surface brightness 
of the disk /x<j (0) as a function of ra- A correlation over- 
all is seen, with fainter central brightness corresponding 
to a larger disk scale length, but there is a clear offset 
between disks of pseudobulge hosts and those of classical 
bulge hosts. For a given r<j, the disk hosting a pseudob- 
ulge is fainter than the one hosting a classical bulge. 

We have performed two-sample Kolmogorov - Smirnov 
tests to compare the distributions of central surface 
brightness, absolute magnitude of disk, and r& for clas- 
sical and pseudobulge hosting SOs. We find that these 
samples of classical bulge hosts and pseudobulge hosts 
could not have been drawn from the same parent popu- 
lation, with at least 99.9% confidence. 

In his paper, Gadotti (2009) finds that the disks of 
pseudobulge hosts are more extended and have a fainter 
surface brightness compared to those of classical bulge 



hosts but the overlap of the two kinds of galaxies is sig- 
nificant. It is important to note that Gadotti (2009) fo- 
cusses on bulge dichotomy for a sample comprising of a 
mixture of different morphological types while our sam- 
ple comprises SOs alone. 

4. DISCUSSION 

We have presented the first systematic study of pseu- 
dobulges in SO galaxies with emphasis on signatures of 
their evolutionary processes on their progenitor disks. 
We use the position of the bulge on the Kormendy dia- 
gram as an initial classification criterion for determining 
the nature of the bulge. To make our classifications more 
secure we have also used the criterion proposed by Fisher 
& Drory (2008) which involves using a division line of 
n = 2 on the Sersic index scale. There are 27 pseudob- 
ulges in our sample of which two belong to bright SOs. 
Thus pseudobulges occur preferentially in fainter galaxies 
(Mr > —22.66, AB system). Using plots between bulge 
effective radius versus disk scale length and disk central 
surface brightness versus disk scale length, we demon- 
strate that distinct trends are followed by classical and 
pseudobulge host SO galaxies. 

The disks of pseudobulge hosting SOs have on average 
smaller scale lengths, lower central surface brightness and 
luminosity. One can either interpret this as pseudobulges 
preferentially occurring in a different population of disks 
or that these signatures are an imprint of the processes 
that drove the pseudobulge growth. It is not clear how to 
explain this finding using one of the known mechanisms 
for pseudobulge formation. For example, transforming 
a spiral with a pseudobulge via ram pressure stripping 
(Abadi et al. 1999) can perhaps cause disk fading but 
may not cause a change in the scale length resulting in 
the lower rd reported here. Perhaps a combination of 
processes may be able to explain the findings reported in 
this Letter. 

In a future work, we will present a study of the current 
sample in greater detail and search for further signatures 
of the underlying processes involved in bulge and pseu- 
dobulge formation in SO galaxies. 
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